Dual chiral density wave in quark matter 
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Abstract 

Possible manifestation of a dual chiral density wave in quark matter is discussed in relation to 
the chiral symmetry of QCD, which is described by a dual standing wave of the scalar and pseudo- 
scalar densities. It is demonstrated that quark matter is unstable for forming the dual chiral density 
wave at moderate densities in the Nambu-Jona-Lasinio model; accordingly the critical density for 
restoration of chiral symmetry becomes higher than ever considered. A unique magnetic property 
is also pointed out in the new phase. 
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Recently condensed matter physics of QCD has been an exciting area in nuclear physics. 
Superconductor model of the vacuum is a classic one, firstly considered by Nambu and 
Jona-Lasinio to understand the nucleon mass Nowadays it is believed that quark-anti- 
quark pairs condense in the vacuum, leading to the spontaneous symmetry breaking (SSB) 
of chiral symmetry in QCD. In the past decade the color superconductivity (quark-quark 
pair condensation) in quark matter has been extensively studied by many authors [2] . There 
are also some works to inquire the coexistence of these two phases on the phase diagram of 
QCD 0. 

On the other hand, the possibility of ferromagnetism (FM) in quark matter has been 
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also discussed jj, |5|| , stimulated by the discovery of magnetars p . It has been suggested 
^ that ferromagnetic phase exists at rather low densities in analogy with an electron gas 
discussed by Bloch fl] and demonstrated that relativistic effects play an important role in 
quark matter. We have also studied a coexistent mechanism of FM and CSC by assuming 
the spin-parallel pairing within the mean-filed approximation j^. 

It should be noted that CSC arises from particle-particle (p-p) correlation, while particle- 
hole (p-h) correlations should be responsible for SSB of chiral symmetry and FM. Besides, 
different types of the p-h condensations has been proposed 0, Q] in relation to chiral sym- 
metry, where the scalar or tensor density forms a standing wave with a finite momentum, 
called the chiral density wave. The instability for the density wave in quark matter was first 
discussed by Deryagin et al. jnl at asymptotically high densities where the interaction is 
very weak, and they concluded that the density-wave instability prevails over the BCS one 
in the large Nc (the number of colors) limit due to the dynamical suppression of colored 
pairing interactions. 

We consider here another type of the density wave due to chiral symmetry in QCD. 
In this letter we are concentrated on the flavor-Sf/ (2) case, and study the flavor- singlet 
scalar density, (gg), and the flavor-triplet pseudo-scalar density, {qi'j^Tiq), which transform 
as the [1/2,1/2] representation in SU{2)l x SU{2)r chiral symmetry. In the SSB phase, 
they satisfy the constraint, (gg)^ + (g^75Tjg)^ = in the chiral limit within the mean- 
field approximation, and the ground states are degenerate on the hypersphere in the chiral 
space, prescribed by the scalar and pseudo-scalar mean-fields, if both densities are spatial 
constants. 

We introduce a dual chiral- density wave (DCDW) state, where both the scalar and 



2 



pseudo-scalar densities always reside on the hypersphere with a constant modulus A, while 
each density is spatially non-uniform; restricting ourselves to the electric charge eigenstate, 
we consider the following configuration in quark matter, 



(qq) = Acos6'(r) 
{qilbT^q) = Asm9{r). (1) 

The chiral angle 6 represents the degree of freedom of the Nambu-Goldstone mode or the 
directional mode on the hypersphere. Taking the simplest but nontrivial form for the chiral 
angle 6 such that ^^(r) = q • r, we call this configuration DCDW. 

It should be noted that we can construct the DCDW state by acting a space-dependent 
chiral transformation such as q — > exp(z75^(r)/2)g, on the usual SSB phase where only the 
scalar density condenses. This is nothing but a kind of Weinberg transformation jl^ . 

When the chiral angle has some spatial dependence, there should appear two extra terms 
in the effective potential as consequences of the local chiral transformation: one is the 
interaction term, ^(757 ■ V6/2)q, due to the non-commutability of 6{r) with the kinetic 
(differential) operator in the Dirac operator. It may be easily seen that V6 acts as an 
"external" axial-vector field for quarks and resolves the degeneracy of the single-particle 
energy spectra for different spin states. In this case, by a suitable rearrangement of the two 
Fermi seas we can always find an optimal configuration to lower the total energy. Another 
one is nontrivial and comes from the vacuum-polarization effect: the energy spectrum of 
quarks is modified in the presence of ^^(r) and thereby the vacuum polarization should give 
an additional term, oc (V6')^ in the lowest order. This can be regarded as an appearance of 
the kinetic energy term for DCDW through the vacuum polarization 3|- Thus, the Fermi 
sea works for DCDW, while the Dirac sea against it; when the interaction energy is superior 
to the kinetic energy, quark matter becomes unstable to form DCDW. 

We will see that the mechanism is quite similar to that for the spin density wave suggested 
by Overhauser 16, ^|, and entirely refiects many-body effects. In general, such a density 
wave is favored in 1-D (one spatial dimension) systems with the wave number |q| = 2pF 
{Pp : Fermi momentum) according to the Peierls instability j^, e.g., charge density waves 
in quasi- 1-D metals The essence of its mechanism is the nesting of Fermi surfaces and 
the level crossing of the single-particle energy spectra with a relative momentum q. In the 
higher dimensional systems, however, the nesting is incomplete and the density wave should 



be formed provided the interaction of a relevant (p-h) channel is strong enough. For the 3-D 
electron gas, it was shown by Overhauser [ly, |l7| that the paramagnetic state is unstable 
with respect to the formation of the static spin-density wave, in which the energy spectra 
of up- and down-spin states with a finite relative momentum exhibit a level crossing and it 
is resolved by the spin exchange interaction, while the wave number does not coincide with 
2pF completely because of the higher dimension. 

We explicitly demonstrate that quark matter becomes unstable for a formation of DCDW 
at moderate densities, using the Nambu-Jona-Lasinio (NJL) model which is originally one of 
the schematic models to describe realization of chiral symmetry in the vacuum []| . Recently 
the model has been also considered as an effective model of QCD embodying SSB of chiral 
symmetry in terms of quark degree of freedom [isj ]. 

We start with the NJL Lagrangian with Nf = 2 flavors and Nc = 3 colors, 

Cnjl = - m,)^ + G[{iji}f + (t^ztst^)'], (2) 

where rric is the current mass, rrtc — 5MeV. Under the mean-fleld approximation, we keep 
the pseudo-scalar as well as the scalar mean-fleld. In the usual treatment to generate a 
dynamical quark mass and see the restoration of chiral symmetry within the NJL model, 
the pseudo-scalar mean-fleld is implicitly discarded. However, note that this is justifled only 
for the vacuum since it is the deflnite eigenstate of parity, and there is no compelling reason 
at flnite baryon density. We assume here the following form for the mean-flelds, 

{ipip) = A cos(q ■ r) 
(7/^275X3?/') = Asin(q-r), (3) 

in the direct channel. Accordingly, we deflne a new quark fleld ipw by the Weinberg trans- 
formation, 

= exp[ii75r3q ■ r/2]^, (4) 

to separate the degrees of freedom of the amplitude and the phase of DCDW in the La- 
grangian. In terms of the new fleld the effective Lagrangian renders 

UiF = ^w[i9 -M- 1/275X3^]^^/ - + £sB, (5) 
where £sb is a small residual term due to the current mass mc, 

£sB = —mcV'iy [exp(i75r3q ■ r) — Ij^/'iy, and we put M = rric — 2G/S. and = (0, q). In the 



following, we take the chiral limit {rric = 0) discarding £sb- The form given in (jSJ appears 
to be the same as that in the superconductor model of the vacuum except an "external" 



axial- vector field, q, generated by the wave vector of DCDW [18| ; the amplitude of DCDW 
generates the dynamical quark mass in this case. Accordingly we can see that our theory 
becomes trivial in the chiral-symmetry restored phase. 

The Dirac equation for ipw then gives a spatially uniform solution, ipw = 
uw{p) exp(zp ■ r), with the eigenvalues 



E^ip) = \jEl + |q|V4 ± v^(p . q)2 + M^jq^, E, = {M' + jpH^/^ (6) 

for positive energy (valence) quarks with different polarizations denoted by the sign ±. For 
negative energy quarks in the Dirac sea, they have an energy spectrum symmetric with 
respect to null line because of charge conjugation symmetry in the Lagrangian (0). 

The single-particle energy spectra (jH)) exhibit an analogue of the exchange splitting be- 
tween energies with different spin states in the Stoner model Hereafter, we choose q/ /z, 
q = (0, 0, q) with g > 0, without loss of generality. Note that we need not distinguish two 
flavors in the energy spectra, while their eigenspinors are different to each other. The energy 
spectra © show a salient feature: they break rotation symmetry due to the couping term 
of the momentum and the wave vector, and lead to the axially-symmetric deformed Fermi 
seas. 

The effective potential then reads 



= fival + f^vac + My AG (7) 



aotai = 7/7^ [{e; - + {e; - ^,)e^] -jJ^^[e; + e;] + (m - m^y/AG 



after summing up all the energy levels, where 6± = d{fi — E^), fi is the chemical potential 
and 7 the degeneracy factor 7 = NfNc{= 6). The first term fivai is the contribution by the 
valence quarks filled up to the chemical potential fi in each Fermi sea, while the second term 
fivac is the vacuum contribution that is formally divergent. We shall see both contributions 
are indispensable in our context. Since the NJL model is nonrenormalizable, we need some 
regularization procedure to extract a meaningful finite value for the vacuum contribution. 
Note that we cannot apply the usual momentum cut-off regularization (MCOR) scheme to 
f2vac, since the energy spectrum has no more rotation symmetry. Instead, we adopt the 
proper-time regularization (PTR) scheme which may be a most suitable one for our 



purpose, since f^vac counts the vacuum-polarization contributions under the "external" axial- 
vector field q. Incidentally the vacuum polarization effect under the external electromagnetic 
field has been treated in a gauge invariant way by using the PTR scheme, where the energy 
spectrum is also deformed depending on the field strength [iQ^ . 
The vacuum contribution fivac can be represented as 

/d'^p 
j^--^tT\n{Sw{p)/So{p)), (8) 

where Sw{p) is the quark propagator, Swip) = ~ M — 1/275X3^]"^. Here we subtracted 
an irrelevant constant Qrei = il ! 'iV/(2vr)'^trln(5'o(p)) with Sq{j)) = — rn^dY^ to make 
the following procedure mathematically well-defined, which is the quark propagator with an 
arbitrary mass mref . Introducing the proper-time variable r, we eventually find 



^ 'J f°° dr f°° dp 



87r3/2 7g ^5/2 J__^ 2tx 



aef, (9) 



which is obviously reduced to the standard formula [13] in the limit g — > (no DCDW). 
The integral with respect to r is not well defined yet as it is, since it is still divergent due 
to the r ~ contribution. Regularization further proceeds by replacing the lower bound of 
the integration range by l/A^, which corresponds to the momentum cut-off in the MCOR 
scheme. 

Now we examine a possible instability of quark matter with respect to formation of 
DCDW. For the vacuum contribution, we can easily see fi^ac > ^vac|q-+0) which means the 
vacuum (the Dirac sea) is stable against formation of DCDW, as it should be. For small g, 

f^vac = f^vac|g-.o + J—^J{M^|^^)q^ + 0{q') (10) 

where J{x) is a function, J{x) = — a;Ei(— a;), with the exponential integral Ei(— x). The 
coefficient of term in Q^ac, which is called the spin stiffness parameter /3vac, has a definite 
physical meaning: since the pion decay constant is given in terms of J{x) within the NJL 
model Q, 
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= ^JiMVli% (11) 
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/3vac can be written as /3vac = ^fw- There are two remarks in order. One is that there appears 
no divergence as A ^ cxo in the higher order terms in q^. Hence, strictly speaking, higher 
order terms should be discarded [l3|, since the divergent term should dominate over other 
finite terms. However, we keep them here to preserve self-consistency after introduction 



of the regularization in the nonrenormahzable theory. The other is that Eq. (jlUp suggests 
that the vacuum-polarization effect precisely provides the expected kinetic-energy term for 
DCDW, except a "running" pion decay constant (fTTj) . 

For given /i, M and q we can evaluate the valence contribution fi^ai using Eq. (jUj) and 
write down the general formula analytically. For small q we have a following expansion, 

^^vai = ^^vaiUo - ^M^q'H{^/M) + 0{q') (12) 



in 



where fivail^^o = 7/(247r2) /xV/^^ - (SM^ - 2/i2) - 3M^ln(/i + V/i^ - MVm) 



for nor- 



mal quark matter and H{x) = ln{x + va?— T). The spin stiffness parameter then reads 

/5vai = -£^M'H{fi/M) (13) 

for valence quarks. Since the function H{x) is always positive and a monotonously increasing 
function, we can see /3vai < 0. 

The total spin-stiffness parameter, Ptot = Pvai + Pva.c, depends on the dynamical mass and 
the chemical potential for a given A. In the vacuum it is positive, and the necessary condition 
to form DCDW then reads f3tot < due to the finite-density effect. When this condition 
is fulfilled, the optimal value of q is determined by the minimum of the effective potential 
((Tj). We can see that the optimal value of q is 0{2pp) with the quark Fermi momentum 
Pp without recourse to explicit calculations. First, consider the energy spectra for massless 
quarks (see Fig.l). As is already discussed, our theory becomes trivial in this case and we 
find two energy spectra 



E^{p) = ^pl + {\p,\±q/2y, P± = (P..,P„0), (14) 

essentially equivalent to the usual ones E'^i^p) = |p|(= E{p)); taking a linear combination 
of the eigenspinors ^lJ^, we can reconstruct the eigenstates of definite chirality 75 to get the 



energy spectra, E'^ip) = yp^ + {pz ± q/'^Y- Then we can see a level crossing on the = 
plane. Once the mass term is taken into account this degeneracy is resolved and the energy 
gap appears there. The resultant spectra are given by ©. Hence we have always an energy 
gain by filling only the lower-energy spectrum E~ (p) up to the Fermi energy, if the relation 
,g = 0{2pf), holds. Thus, we find that this mechanism is very similar to that of spin density 
wave by Overhauser Q . 

Taking the extremum of the effective potential (|7j) with respect to the order-parameters M 
and q, we can determine their values for given baryon-number densities. Fig. |21 demonstrates 



-q/2 q/2 

Pz 

FIG. 1: Energy spectra for pj^ = 0. with M = (thick sohd and dashed hnes). with the 
definite chirahty is also shown for comparison (dotted hne). We can see there is a degeneracy of 
E^ at ^2 = for M = 0, while it is resolved by the mass (thin solid and dashed lines). 

the behaviors of the optimum values as functions of /i at zero temperature for a parameter 
set, GA? = 6, A = 860 MeV. It is found that the magnitude of q becomes finite just before 

0.8 
0.6 
0.4 
0.2 
0.0 

0.46 0.48 ^ol 0.50 0.52 *^c2 0.54 

FIG. 2: The wave number q and the dynamical mass M arc plotted as functions of the chemical 
potential at zero temperature. Solid (dotted) line for M with (without) DGDW, and dashed line 
for q. 

the ordinary chiral-symmetry restoration, and DCDW survives in the finite range of 
A^ci < A* < A*c2, which corresponds to the baryon-number densities Ph/ Po — 3.62 — 5.30. 
The dynamical mass remains finite in this density region, which means chiral restoration is 
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delayed in the presence of DCDW. The wave number q increases with /i, while its value is 
smaller than 2pf{— 2/i for massless quarks) due to the higher dimensional effects. Actually, 
the ratio of the wave number and the Fermi momentum (at normal phase g = M = 0) lies 
in the range, q/pr = 1.17 — 1.47. We thus conclude that DCDW is induced by finite-density 
contributions, and has the effect to extend the SSB phase (M ^ 0) to high density region, 
which suggests another path for chiral-symmetry restoration by way of the DCDW state at 
finite density. 

Here we would like to point out an interesting aspect of the DCDW state. The quark- 
DCDW coupling is spin dependent and we can see that it gives rise to a spatial oscillation 
of the magnetic moment in quark matter. The magnetic moment is given as 

M, oc (^cTisV') = Mf cos(q ■ r), (15) 

where MY {= {'4^w(^i2''Pw)) is proportional to the dynamical mass. Thus DCDW can be 
regarded as a kind of the spin density wave jlfil . u\ and should have phenomenological 
implications on the magnetic properties of compact stars. 

In this Letter we have only considered the direct channel (Hartree term) of the four 
Fermi interaction. When we consider the exchange contributions (Fock term), there appear 
additional quark-quark interactions in the vector and axial-vector channels by way of the 
Fierz transformation . It is well known that the vector contribution is rather trivial 

and summed up as a renormalization of the chemical potential. We are then interested in the 
axial- vector contribution in our context. We may see that it renormalizes the wave vector 
of DCDW; this can be regarded as a kind of vertex renormalization of the quark-DCDW 
coupling. This subject will be discussed elsewhere 

It would be interesting to recall that DCDW is similar to pion condensation within the 
a model, considered by Dautry and Nyman [2^, where a and vr*^ meson condensates take 
the same form as Eq. ([T)). So it might be possible to connect pion condensation before 
deconfinement with DCDW after it by a symmetry consideration. However, the restoration 
of chiral symmetry, especially the vacuum contribution, has been poorly considered in the 
context of pion condensation. 

Finally, it might be interesting to generalize the DCDW configuration (jHl) by incorporat- 
ing the radial degree of freedom; one may generalize the amphtude to be spatially dependent 
as well a la refs.]^, but one must consider a genuine non-uniform matter in this case. 
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Anyhow, if two degrees of freedom work coherently, we may have more favorable configura- 
tions than considered here. 

We thank T. Maruyama, H. Yabu and K. Nawa for discussions. This work is supported 
by the Grant-in-Aid for the 21st Century COE "Center for the Diversity and Universality 
in Physics " from the Ministry of Education, Culture, Sports, Science and Technology of 
Japan. It is also supported by the Japanese Grant-in-Aid for Scientific Research Fund of 
the Ministry of Education, Culture, Sports, Science and Technology (13640282, 16540246). 



[1] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345. 
[2] D. Bailin and A. Love, Phys. Reports 107(1984) 325. 

For a review, e.g., M. Alford, Ann. Rev. Nucl. Part. Sci. 51 (2001) 131. 
[3] J. Berges and K. Rajagopal, Nucl. Phys. B538 (1999) 215. 

R. Rapp, T. Schaefer, E.V. Shuryak, M. Velkovsky, Phys. Rev. Lett. 81 (1998) 53. 
[4] T. Tatsumi, Phys. Lett. B489 (2000) 280. 



[5] A. Niegawa, |hep-ph/ 0404252( 



[6] For a recent review, P.M. Woods and C. Thompson, astro-ph/0406133 
[7] F. Bloch, Z. Phys. 57 (1929) 545. 

K. Yosida, Theory of magnetism (Springer, 1998). 
[8] E. Nakano, T. Maruyama and T. Tatsumi, Phys. Rev. D68 (2003) 105001. 
[9] E. Shuster and D. T. Son, Nucl. Phys. B573 (2000) 434. 
[10] B.-Y. Park, M.Rho, A.Wirzba and I.Zahed, Phys. Rev. D62 (2000) 034015. 

R. Rapp, E.Shuryak and I. Zahed, Phys. Rev. D63 (2001) 034008. 
[11] D.V. Deryagin, D. Yu. Grigoriev and V.A. Rubakov, Int. J. Mod. Phys. A7 (1992) 659. 
[12] S. Weinberg, The quantum theory of field //(Cambridge, 1996). 
[13] T. Eguchi and H. Sugawara, Phys. Rev. DIO (1974) 4257. 

K. Kikkawa, Prog. Theor. Phys. 56 (1974) 947. 
[14] R. E. Peierls, Quantum Theory of Solids (Oxford University Press, London, 1955). 
[15] G. Griiner, Rev. Mod. Phys. 60 (1988) 4. 

[16] A.W. Overhauser, Phys. Rev. Lett.4 (1960) 462; Phys. Rev. 128 (1962) 1437. 
[17] G. Griiner, Rev. Mod. Phys. 66 (1994) 1. 



10 



[18] S.P. Klevansky, Rev. Mod. Phys. 64 (1992) 649. 

T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994) 221. 
[19] J. Schwinger, Phys. Rev. 92 (1951) 664. 

[20] T. Maruyama and T. Tatsumi, Nucl. Phys. A693 (2001) 710. 

[21] E. Nakano and T. Tatsumi, in preparation. 

[22] F. Dautry and E.M. Nyman, Nucl. Phys. 319 (1979) 323. 

K. Takahashi and T. Tatsumi, Phys. Rev. C63 (2000) 015205; Prog. Theor. Phys. 105 (2001) 

437. 

M. Kutschera, W. Broniowski and A. Kotlorz, Nucl. Phys. A516 (1990) 566. 
[23] It has been shown that similar mechanism by valence quarks also works in another context of 
the chiral density wave P, Q] . 



11 



